Superfluidity of trapped dipolar Fermi gases. 
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We derive the phase diagram for ultracold trapped dipolar Fermi gases. Below the critical value 
of the dipole-dipole interaction energy, the BCS transition into a superfluid phase ceases to exist. 
The critical dipole strength is obtained as a function of the trap aspect ratio. The order parameter 
exhibits a novel behavior at the criticality. 
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The quest for the superfluid (BCS) transition J,J in 
trapped Fermi gases is one of the most challenging goals 
of modern atomic, molecular and optical physics. The 
possibility of the BCS transition for trapped gases with 
attractive short range interactions has been predicted in 
and has been a subject of extensive experimental re- 
search since then. In typical experiments evaporative 
cooling is used to cool fermions. However, as the Pauli 
principle forbids the s-wave scattering for fermions in the 
same internal state, Fermi-Fermi jj| or Fermi-Bose 
mixtures have to be used to assure coUisional thermaliza- 
tion of the sample. Such a combination of evaporation 
and sympathetic cooling allows to reach temperatures 
T ~ O.lTp, where Tp is the Fermi temperature at which 
the gas exhibits quantum degeneracy. Unfortunately, 
predicted critical temperatures for the BCS transition, 
Tc, are much smaller than Tp. One way to circumstance 
this difficulty is to increase Tc] this way may be achieved 
by increasing the strength of the atomic interactions em- 
ploying a Feshbach resonance, which allows to make the 
atomic scattering length Os large negative. Such "reso- 
nance superfluidity" should lead to superfluid transition 
at Tc ~ O.lTp 5]. Another way is to use the cooling 
scheme that can overcome the effects of Pauli blocking, 
such as appropriately designed laser cooling j^]. Yet an- 
other promising route is to locate the ultracold gas in an 
optical trap and enter the strongly correlated regime by 
controlling the lattice potential 7\ . 

The temperature of the BCS transition in a two- 
component Fermi gas, however, depends also drastically 
on the difference of the concentrations of the two com- 
ponents, which presents another experimental obstacle. 
This problem is not relevant for a polarized Fermi gas 
with long-range interactions, such as dipole-dipole inter- 
actions, and that is why there has been a considerable 
interest recently in studying the BCS transition in dipo- 
lar Fermi gases. The possibility of the Cooper pairing 
has been predicted in Refs.0. The critical temperature 
(including many-body corrections) and the other param- 
eters have been obtained in Ref.0- At this point it 
is worth mentioning that possible realizations of dipo- 
lar gases include ultracold heteronuclear molecular gases 
[inj l , atomic gases in a strong DC electric field , atomic 



gases with laser- induced dipoles or with magnetic 
dipoles 01 . For dipolar moments d of the order of one 
Debye and densities n of lO^^cm"'^, Tc should be in the 
range of lOOnK, i.e. experimentally feasible. 

Dipole-dipole interaction is not only of long-range, but 
also anisotropic, i.e. partially attractive and partially 
repulsive. Thus for trapped gases, the nature of the in- 
teraction may be controlled by the geometry of the trap. 
For a dipolar Bose gas in a cylindrical trap with the axial 
(radial) frequency LOz^cOp), there exist a critical aspect ra- 
tio A = {ujz/ujpy^'^ , above which the Bose-condensed gas 
collapses if the atom number is too large |Q| , and below 
which the condensate exhibits the roton-maxon instabil- 
ity 01 . The effects of trap geometry are expected to 
dominate also the physics of trapped dipolar Fermi gases 
(16||. So far, we have only been able to derive analytic cor- 
rections to Tc in "loose" traps, and to consider the case 
of an infinite "slab" {ojp = 0, is finite). In this case 
there exists a critical frequency above which the dipole 
interaction is predominantly repulsive, and the superfluid 
phase ceases to exist. 

In this Letter we present the ultimate solution of the 
fundamental problem of the effect of trap geometry on 
the BCS transition in trapped dipolar Fermi gases. We 
calculate the phase diagram in the plane F — A~^, where 
F ~ nd? / IX is the dipole-dipole interaction energy per 
particle in the units of the chemical potential /i. Be- 
low the critical value F < Fc, the BCS transition does 
not take place. We determine dependence of Fc on A, 
and calculate the order parameter at the criticality. The 
order parameter exhibits a novel oscillatory behavior in 
strongly elongated cylindrical traps. 

We consider a dipolar single-component Fermi gas in 
a cylindrically symmetric trap described by the Hamilto- 
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where m is the atoms mass, T4rap(r) = m[ujp{x^ -f y^) + 
w^^^J the trapping potential, /i the chemical potential. 
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and Vdip(r) = {d'^/r^){l - Sz'^/r'^). The dipoles are 
assumed to be polarized along the z-direction, ^p^r) 
and ipir) are the atomic creation and polarization op- 
erators that fulfill canonical anticommutation relations. 
Our aim is to apply the BCS theory (see. e.g. 0) 
to the system described by Eq. 0J, calculate the crit- 
ical temperature, and the superfluid order parameter 

A(r, r') = Vdip{r — r') (^l){r)■^l){r')^. In particular, we 

are interested in the critical value of the aspect ratio A 
below which the BCS pairing does not take place for a 
given strength of the dipole interaction. The BCS the- 
ory leads to the gap equation which is a linear integral 
equation for A(r, r') (a function of 6 variables) with a 
kernel that depends on 12 coordinates. Numerical solu- 
tion of this equation is impossible without further an- 
alytical simplifications. As shown in Ref.j^, the domi- 
nant contribution to the BCS pairing comes from the p- 
wave scattering (contributions of higher angular momen- 
tum, although present due to the long-range character 
of the dipole-dipole interaction, are numerically small) 
with zero projection of the angular momentum on the 
X-axis, Iz = 0, in which the interaction is attractive (see 
also Ref. 0). In the p-wave channels with = ±1 
the interaction is repulsive and, therefore, results only in 
renormalizations of a Fermi-liquid type. The latter are 
controlled by the small parameter of the theory F and 
will be neglected. Thuerefore, for the considered pairing 
problem we can model the dipole-dipole interaction by 
the following (off-shell) amplitude 
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where p is the incoming momentum, p' the outgoing one, 
and 7i(-E) some function of the energy E. The ampli- 
tude Fidescribes anisotropic scattering only in the p- 
wave channel with the projection of the angular momen- 
tum Iz = 0. The function 7i(£') obeys the equation 
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that follows from the Lipmann-Schwinger equation for 
the off-shell scattering amplitude, and is assumed to be 
negative, ji{E) < 0, in order to guarantee the BCS pair- 
ing. The cut-off parameter A ensures the convergence of 
the integral and, in fact, can be expressed in terms of the 
observable scattering data corresponding to on-shell scat- 
tering with p — p' and E ~ p^ /m. Tt follows from the 
above equation that 71 (-E) is inversely proportional to E^ 
71 (i?) = 71 {2mE)~^, with some coefficient 71. There- 
fore, the on-shell amplitude is energy independent, as it 
should be for low-energy scattering on the dipole-dipole 
potential (see Ref. 17]). 



The coefficient 71 determines the value of the critical 
temperature of the BCS transition in a spatially homo- 
geneous gas and hence, using the results of Ref. 0, can be 
expressed through the dipole moment d. In this case the 
order parameter has the form A(p) — Pz^o with some 
constant Aq, and the linearized gap equation results in 
the equation for the critical temperature Tc 
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where = p^ /2ra — /i and the bare interaction is 
renormalized in terms of the scattering amplitude with 
71 (/i) = ^i/p]? at the Fermi energy = M = p'p/^m 
along the lines of Ref. (pp is the Fermi momentum) . 
After integrating over p we obtain the following equation 
for Tr: 
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where lyp — mpi?/27r^ is the density of states at the Fermi 
energy and C — 0.5772 the Euler constant, and therefore 
71 = — 24(i^/7r in order to reproduce the result of Ref.0 
for T,. 

In the ordinary space, the scattering amplitude Fiis 
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where r and r' are the relative distances between the two 
incoming and two outgoing particles, respectively. There- 
fore, the order parameter in a trapped gas, A(ri,r2) ^ 
(V'(ri)'(/'(r2)), has the form A(ri,r2) = 9z(5(r)Ao(R), 
where r = ri — r2 and R = (ri -f r2)/2, and the cor- 
responding equation for the critical temperature is 
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where ii C(ni), n = (?^x, 't-z) are the 

harmonic oscillator quantum numbers, '^(n) = 
h\ujz{nz ^ ijjp{nx + Uy + \)\ — ^, and 

M„,„,(R) EE Mitn,Az)Mi''^.nUx)Mif,ln,,{y) with 

( 1 

M,\/n^{z) = ^[Lpni{z)dz'Pn2{z) - 'Pn2{z)dzVnAz)], 

Mn^uiix) = ipni{x)'fn2{x), and similar for Aflil^a (y) , 
are the harmonic oscillator wave functions. 

The gap equation Q is still hardly tractable numeri- 
cally. We employ thus the series of additional approxima- 
tions. We assume a large number of particles in the trap 
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and hence, the chemical potential fj, is much larger than 
the trap frequencies, /x ^ Wz,Wp. Then, for the most im- 
portant for the BCS pairing states near the Fermi energy 
ep = fi we can use the WKB approximation for the wave 
functions ipn while calculating the functions Mnin2(R)- 
Another simplification is due to the fact that the BCS 
order parameter Ao(R) varies slowly on an interparti- 
cle distance scale n~-^/^ ~ fi/PF, where pp ~ ^/2mJI is 
the Fermi momentum in the center of the trap in the 
Thomas- Fermi approximation. As a result, the pairing 
correlations are pronounced only between close in energy 
states. Therefore, one can neglect /Am in the denomi- 
nator of the second term in Eq.©, as well as rapidly os- 
cillating contributions in the functions M^j^n^ (R). Under 
these assumptions these functions become proportional 
to the Chebyshev polynomials \Jn^iiz/lzj\j), Tn{x/lpN) 

^j2Nh/nu}r 



and T„{y/lpN), where = ^/2Nh/nuJi = Iq^V^N. 

The critical aspect ratio Ac corresponds to vanishing 
critical temperature. We therefore calculate the gap 
equation in the limit T <C Wi. After very tedious, 
but fully analytical calculations with the substitution 
Ao(R) ^ Ao(r) = Aa{zR^^l,xR'^l,yR'^l), where x, y, 



and denote 



z are now dimensionless |a;|, \y\, \z\ < 1 
the radius of the gas cloud in the i-direction, calculated 
in the Thomas-Fermi approximation, R^rpp — pp/moji, 
we arrive at the gap equation 



^(l-r2)Ao(r) = (l-r2)-V2|ln 
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where T ^ \^i\vf, ^j.{y) ^ ^ - Vtrap(r), and 
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with bn = 2 for n > 0, i5o = 1, and Mg = max(s^,s'^) 
for s — x,y, z. The above equation can be viewed as the 
equation for an extrcmum of the quadratic form 



F[Ao] = ^ Ao(r)[i(r)<5(r - r') - K{r, r')]Ao(r'), 



where L{r) = 3(1 - r'^)/T - (1 - r'^)^^^\n{[2n{r)/ujz] - 
2(4 — ln4)/3}. The extremum can be found numerically 
using the ansatz 
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with rriz, nip > 0. The form F becomes now a quadratic 
form of the unknown coefficients Cm^mpi and the ex- 
tremum corresponds to the eigenvector of the matrix 
Mm^,mp,n^,np of the quadratic form with zero eigenvalue. 
Such an eigenvalue exists only if the interaction F and 
the trap frequencies coi obey certain constraint. The pa- 
rameter F can be written as F = 36n(0)d^/7r/i, where 
n{0) = {2mii)^^^ /Gn^h^ is the gas density in the center 
of the trap, and hence, for a given dipole moment d, de- 
pends only on the chemical potential /i. On the other 
hand, the chemical potential fi and the total number of 
particles A*" in the trap are related as 3A^ = /i'^/w^w^. 
Therefore, after fixing F and M, the product of the trap 
frequencies ujzt^p is also fixed, and the only free parame- 
ter left is the trap aspect ratio A, which critical value Ac 
could be determined from the above constraint. 

The problem is therefore reduiced to finding the set 
of coefficients Cm^mp and the aspect ratio A such that 
for a given F and N the lowest eigenvalue of the matrix 
Mm^,mp,n^.np of the quadratic form F is zero. The matrix 
M is symmetric, and hence, there exists a solution. 

The calculation of the matrix elements Mm^.mp,n^,np 
is naturally divided into two parts (see Eq. ©). The 
local part with the kernel L{r) is just a 3-dimensional 
integral that can be easily computed using, for instance, 
the Monte Carlo integration routines, such as the VE- 
GAS algorithm from the GSL library [ijj • The non-local 
part with the kernel K{r,r') is a triple sum, which ele- 
ments are 8 dimensional integrals. The straightforward 
approach using the same numerical method as described 
above fails in this case, because it is too time consuming. 
To overcome this problem, we perform the integrations 
over r, r'and C for fixed a, [3 and n's using the same 
method as before, and then use a 2-dimensional spline 
interpolation method to interpolate the integrand for the 
last integrations over a and (3. In this way the time re- 
quired to compute the matrix elements of M for a given 
set of physical parameters, was reduced to about 72 hours 
on a standard workstation. 

The results of the calculations are presented in 2 fig- 
ures. Fig.l shows the desired relation betweem F and the 
aspect ratio A. The two curves correspond to two differ- 
ent numbers of particles. As it could be expected, for 
the larger number of particles, the critical aspect ratio 
Ac is smaller because the interaction is stronger. For a 
pancake trap, A~^ < 1, the interaction is predominantly 
repulsive, and higher values of F for fixed A are required 
to achieve the BCS transition. On the other hand, for 
a cigar trap, A~^ > 1 , the interaction is predominantly 
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FIG. 1: The critical Tc as a function of the aspect ratio A: 
above the depicted curves BCS takes place. The upper (lower) 
curve corresponds Xo N = lO" and N — 2 x 10®. 

attractive and the occurrence of the BCS transition re- 
quires smaller values of F. 

Fig. 2 shows the order parameter Ao(r) for the cigar 
trap with A^^ = 2.2. Amazingly, the order parameter is a 
nonmonotonic function of the distance from the trap cen- 
ter, in contrast to the case of the BCS order parameter in 
a two component Fermi gas with short range interactions 
llfll. This effect persists, although being less pronounced, 
for the case of a pancake geometry. In the axial direction, 
the order parameter A(z, p = 0) develops a minimum at 
p < 1, whereas in the axial direction A(z = 0,p) be- 
comes negative in the outer part of the cloud. This is a 
completely new behavior originated from the anisotropy 
of the interpaticle interaction, that can have profound 
consequences for the form and spectrum of the elemen- 
tary excitations. We expect an appearance of excitations 
with wave functions concentrated in the inner region of 
the atomic cloud, where A ~ 0. This problem, however, 
goes beyond the scope of this paper. 

Summarizing, we have determined the phase diagram 
for trapped dipolar Fermi gases in the F — A~^ plane, 
where F measures the dipole strength and A is the trap 
aspect ratio. The BCS transition at finite temperature T 
is possible iff F > Fc(A). We have calculated the critical 
line Fc(A), and the order parameter at criticality. 
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FIG. 2: Order parameter for the aspect ratio A~^ = 2.2 
(corresponding to cigar geometry). The solid line shows 
Ao(z, p = 0), and the dashed line corresponds to Ao(2 = 0, p). 
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